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A present challenge in complex systems is to identify mesoscale structures that have distinct 
dynamical implications. In this paper we present a detailed investigation of a previously observed 
dynamical equivalence of certian ecological food webs. We show that this equivalence is rooted 
in mesoscale symmetries that exist in these webs. Certain eigenvectors of the Jacobian describing 
dynamical modes of the system, such as specific instabilities or responses to perturbations, localize 
on these symmetric motifs. On the one hand this means that by removing a symmetry from the 
network one obtains a system which has identical dynamics except for the removal of the localized 
mode. This explains the previously observed equivalence. On the other hand it means that we can 
identify dynamical modes that only depend on the symmetric motif. Symmetric structures thus 
provide an example for mesoscale network motifs having distinct and exact implications for the 
dynamics. 
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I. INTRODUCTION 

Making progress in the investigation of complex sys- 
tems requires finding concepts that allow reducing the 
complexity, while leaving emergent-level features of the 
system intact. Networks offer a framework facilitating 
this reduction PHI]. By representing the system as a 
set of discrete nodes connected by links, networks ne- 
glect much of the internal complexity of the system's 
constituents, but retain the complexity of their interac- 
tion. Thereby, they offer a suitable intermediate level of 
description on which a qualitative understanding of the 
systems dynamics can be gained. 

In the present paper we consider specifically the ex- 
ample of ecological food webs Q, i.e. the networks of 
who-eats-who in ecology. In these networks, the nodes 
represent distinct populations, whereas the links repre- 
sent predator-prey interactions. Food webs offer a highly 
simplified description of ecological systems, neglecting for 
instance non-predatory interactions and intra-population 
dynamics. Yet, they capture the complexity of the net- 
work of predatory interactions that forms the backbone 
of most natural ecosystems. 

A detailed analysis of food web dynamics, typically 
comprising 50-500 populations on very different time 
scales, poses a considerable challenge. Therefore, one 
often considers coarse-grained models in which different 
populations are grouped and represented by one aggre- 
gate variable. A classical example are food chain models 
in which whole layers of the food web such as phyto- 
plankton, zooplankton, or fish are described by a single 
variable, respectively [fllzj]- More biological information 
is preserved by aggregating so-called functional groups 
or guilds comprising species which are similar in mor- 
phology, life strategy, or niche 0, S|. Alternatively, it 
has been proposed to aggregate those species that hold 
similar topological positions in the network j9l-[l^. 

From a theoretical perspective, the problem of coarse- 
graining networks is interesting because it touches upon 



mesoscale network properties. We seek to identify a set 
of nodes that are sufficiently similar, such that their ag- 
gregation induces only minor changes in the dynamics. 
The mesoscale properties of food webs have been previ- 
ously studied in terms of community detection [ll| - |l3j 
and motif distributions [T3 - [l6j . Although both commu- 
nities and motifs are known to have a strong impact on 
the dynamics of the system, the precise nature of this 
impact is not known. 

The ultimate goal of coarse-graining is finding a re- 
duced representation of the network that conserves the 
dynamics of the system. This goal is generally hard to 
achieve as it seems to imply knowledge of the dynamics 
of the original system. However, in a recent publication 
[ItJ it was shown in a class of generalized models [l||, [l9| 
that under certain conditions an aggregation is possible 
that conserves the local dynamics around steady states 
and at least certain features of the global dynamics ex- 
actly. Although the previous paper demonstrated this 
aggregation in several examples, it failed to provide a de- 
tailed analysis of the underlying symmetries that makes 
the aggregation possible. 

In the present paper we build on recent advances [2(| 
in the dynamical implications of graph orbits, which cap- 
ture mesoscale symmetries in the network structure and 
extend them to directed and weighted networks. Ap- 
plication to food web dynamics reveals the topological 
mechanism behind the reduction rule identified heuristi- 
cally in Ref. [13 and provides more general rules for food 
web aggregation. We note that the previously observed 
identity is of limited use for the analysis of food web data. 
However, its theoretical investigation provides an exam- 
ple of mesoscale structures having distinct implications 
on the network dynamics. 

The paper starts, in Sec. HH with a brief review of pre- 
vious results. In Sec. IHII we introduce the notion of sym- 
metries in food webs. Furthermore, we show that the 
nodes of a symmetry carry localized dynamical modes 
that characterized solely by the mesoscale structure of 
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their symmetry. In Sec. IIVI we establish the dynami- 
cal modes for several example food web symmetries and 
derive the coarse-graining rules allowing to remove dy- 
namical modes of on the mesoscale. 



II. DYNAMICAL EQUIVALENCE OF FOOD 
WEBS 



The application considered below builds on previous 
results that have been obtained with a generalized model 
(GM) of ecological food webs [HI, [l9j]. In this section 
we provide an overview of the approach of generalized 
modelling and the previous results on food web stability. 

A GM is a system of differential equations in which not 
all processes are restricted to specific functional forms 
[l8j . For instance, the food web model studied below 
describes the dynamics of ./V populations X\ , . . . , Xjy by 
TV differential equation of the form 



dt 



X, = Gi{X) + Si{Xi) - Li(X) - Mi(Xi), (1) 



where Gi, Li, M iy and Si arc unspecified functions de- 
scribing the gain by predation (Gi), the loss by the 
predation(Li), the loss due to mortality (Mi), and the 
gain by reproduction (Si). 

Generalized models are typically too general to com- 
pute the number or location of steady states without 
further assumptions. Yet, they can be analysed by lo- 
cal stability analysis. For this purpose, we assume that 
a steady state X* — (X^ , . . . , X^) T were known and 
then formally linearize the system around this unspecified 
steady state. The coefficients of the linearized system can 
be expressed by the Jacobian matrix, Jij = 
which depends on the derivatives of the unspecified func- 
tions in the model. For instance, in the food web model 
this unspecified information appears in expressions such 

as §^ . This derivative describes an unknown but con- 

stant quantity or, in other words, an unknown scalar pa- 
rameter of the system. Because this parameter is hard 
to interpret, we use a slightly different parametrization, 
which is obtained either by a special normalization pro- 
cedure [HI or directly by the identity 



dGj 
dX, 



G* dlog d 



X* dlogX, 



(2) 



which is true for G*,X* > (a condition that is generally 
met by definition; the special case of X* = is discussed 
in Ref.iK). 

The expression on the right-hand-side of ([2]) is a prod- 
uct of two factors that can be interpreted directly in most 
applications. The first factor is a per-capita rate. Such 
rates have the dimension of inverse time and can be di- 
rectly interpreted as characteristic turnover rates, i.e. the 
inverse of the mean life expectancy of individuals. The 
second factor is a logarithmic derivative. Such deriva- 
tives are also called elasticities and have been proposed 
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FIG. 1. Equivalence of the bifurcation diagrams of the dy- 
namics in food webs. The bifurcation diagrams of different 
food web models are exactly identical, if plotted as a func- 
tion of the generalized parameters, as shown in the three- 
parameter bifurcation diagrams in the top row (reprinted from 
Ref. [l?]). The parameter space in this example is spanned 
by three generalized parameters denoting the nonlinearity of 
mortality m, the sensitivity of predators to prey abundance 
g, and the ratio between the timescales of predator and prey 
dynamics r. Every point in the parameter space corresponds 
to a specific steady state. The surfaces in the diagrams mark 
the location of bifurcation points of Hopf (red, green) and 
saddle-node (blue) bifurcations. Steady states are stable if 
they lie in the front volume of the bifurcation diagram. The 
two diagrams correspond to different food webs, indicated 
symbolically in the corner. In the symbolic representation 
nodes indicate populations and arrows indicate biomass flows 
due to predator-prey interactions. Using this notation the 
dynamical equivalence conjectured in Ref. Qj] is shown sym- 
bolically in the lower row. 



originally in economic theory [221 ] and are also used in 
the context of metabolic control theory (23|. They can 
be estimated well from data and interpreted straightfor- 
wardly. For every power-law, f(x) = Ax p , the logarith- 
mic derivative is = p, independently of A or x. 
Thus, for instance, any linear function has an elasticity 
of one regardless of the slope. For functions that are not 
power-laws the elasticity still provides an intuitive non- 
linear measure of the sensitivity in the steady state. 

From the parametrized Jacobian the dynamics close to 
the steady state under consideration can be computed. 
The steady state is stable if all eigenvalues of the Jaco- 
bian have negative real parts [24J. Stability is lost when 
a change of parameters causes at least one eigenvalue to 
cross the imaginary axis, leading to a local bifurcation. 
Because the Jacobian is a real matrix the eigenvalues 
in question can be either real or form a complex conju- 
gate eigenvalue pair. In the former case a bifurcation of 
saddle-node type occurs in which the number of steady 
states changes, whereas in the latter Hopf bifurca- 

tion occurs which marks the onset of (possibly transient) 
oscillations. 

The mathematical validity of generalized modelling 
was recently confirmed by rigorous proofs 21 1, whereas 
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applicability to actual problems from different fields was 
demonstrated in several recent publications [lH, I25l43lll. 

A GM for ecological food webs [IH was used in Ref. [T3 
to study critical points in parameter space in which the 
dynamics around steady states changes in bifurcations. 
The authors noted that in the context of the GM food 
webs with different topologies and different numbers of 
species may be described by exactly identical local bifur- 
cation diagrams when these are drawn in the generalized 
parameter space. An example of the observed identity is 
reprinted in Fig. [TJ The figure shows three-parameter bi- 
furcation diagrams for two different food webs for which 
the surfaces formed by bifurcation points of codimension- 
one bifurcations are identical. 

The main conclusion of the previous paper [l7| can be 
summarized in the conjecture that two species that a) 
interact with the same set of topological neighbours and 
b) are described by the same generalized parameters, can 
be aggregated into a single variable without changing the 
local bifurcation points. 

We note that the conjectured identity of bifurcation 
diagrams is not trivially true because of the nonlinearity 
of the underlying equations. Indeed, the identity is much 
harder to observe in conventional models and was hence 
not noted before the publication of Ref. [lj]. 

Finally, consider that the identity extends also to local 
bifurcations of higher codimension (found for instance 
at the intersections of bifurcation surfaces), which have 
implications for the global dynamics of the system 0, 
HH- So the proposed aggregation must also conserve at 
least some features of the global dynamics. 



III. MATHEMATICAL TREATMENT OF 
SYMMETRIES 

The equivalence of food webs described above is contin- 
gent on the presence of a specific symmetry in the topol- 
ogy of the food web. In particular, equivalent reduction 
of the network was found to be possible if at least two 
populations in the network interacted with exactly the 
same neighbours. 

In this section we review and extend some previous 
results on the implications of symmetries in networks. 
Furthermore, we introduce the terminology that is used 
below for establishing a formal link between mesoscale 
symmetries and dynamics. We note that our presenta- 
tion follows in large parts Ref. but extends the theory 
presented therein to directed networks. 

We say that a network contains a mesoscale symmetry 
if there is a permutation of nodes that leaves the adja- 
cency relations unchanged. For instance Fig. [2^) shows 
a small part of a larger network in which two nodes can 
be exchanged without changing the adjacency. We note 
that if multiple symmetries exist they can be combined 
into larger symmetries, but for simplicity we assume that 
the symmetries are minimal, i.e. that they cannot be de- 
composed into symmetries acting on fewer nodes. 
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FIG. 2. Three examples of symmetric structures in networks. 
Orbits are marked by the colouring of the corresponding 
nodes. Simple symmetric structures can consist of a single 
orbit (a, red), or multiple orbits (b, red, yellow), and can po- 
tentially contain nested orbits (c, red, yellow). See text for 
details. 



In the following we denote all nodes that participate 
in a given minimal symmetry as a symmetric structure. 
Further, a set of nodes that are mapped onto each other 
under a minimal symmetry operation are said to form an 
orbit. 

In the simplest case a symmetric structure contains 
only one such orbit (Fig. [2^,). However, some symmet- 
ric structures comprise multiple orbits such that multiple 
permutations have to be carried out in parallel on differ- 
ent sets of nodes to preserve adjacency (Fig.[2}j). 

We refer to the number of nodes belonging to an orbit 
as the size of the orbit. The size corresponds to the num- 
ber of identical graphs generated by the corresponding 
symmetry operation unless a symmetric structure con- 
tains nested orbits (Fig. Oi). 

Although the concept of orbits was originally proposed 
for unweighted undirected networks, it can be intuitively 
extended to weighted and directed ones. In this case we 
additionally demand that symmetries conserve link di- 
rection and weight, such that the corresponding permu- 
tation leaves the matrix of link weights of the network 
unchanged. 

In Ref. it is shown that symmetries in undirected 
networks have a distinct impact on the spectrum of the 
adjacency matrix. In particular, if a symmetric struc- 
ture is present then the eigenvalues can be decomposed 
into two classes. The first class contains the so-called 
redundant eigenvalues. These eigenvalues correspond to 
eigenvectors that are localized on the symmetric struc- 
ture. The second class contains all other eigenvalues. 
They correspond to eigenvectors that are constant on 
the orbits of the symmetric structure. These results sug- 
gest that reducing the network by collapsing all nodes 
in each orbit onto each other leaves the spectrum of the 
adjacency invariant, except for removing the redundant 
eigenvalues. 

In the following we show that very similar results can 
be obtained for the Jacobian matrix of the food web 
model. When the redundant eigenvalues have a nega- 
tive real part the resulting reduction leads to a lower- 
dimensional, but dynamically equivalent, system. To es- 
tablish this reduction we consider the Jacobian as the 
adjacency matrix prescribing the weights and directions 
of the links in a (hypothetical) network. 
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We note that for food webs the network prescribed by 
the Jacobian matrix is not identical to the network of 
predator-prey interactions. For instance, additional mu- 
tualistic interactions may appear between different prey 
populations of a given predator [l8j . However, if different 
populations belong to an orbit in the network of predator- 
prey interactions, and if their generalized parameters are 
identical, then they also form an orbit in the network de- 
fined by the Jacobian matrix. We therefore only have to 
show that the results of Ref. extend to such weighted 
and directed networks. 

Following Ref. [13 we argue that each symmetry gener- 
ates an invariant subspace of the Jacobian matrix. No- 
tably, this subspace contains only vectors that are zero 
outside the symmetry and thus acts as the foundation for 
the set of eigenvectors that are localized on the symmet- 
ric structure. 

We consider a subspace V spanned by vectors v that 
are zero except in the elements belonging to the symmet- 
ric structure. The defining feature of the independent 
subspace is that Jv must lie within V. This condition 
is met when v obeys the "zero-sum" condition for each 
orbit of the symmetric structure, i.e. ^ Vi = 0, where i 
runs over the elements of the orbit. 

To test the localization we first consider a node k out- 
side the symmetric structure. If this node is connected 
to any node within an orbit then symmetry implies that 
it must be connected to all nodes within the same orbit. 
Further, all of these connections must carry the same 
weight (say, 6). We can therefore write the terms of the 
fc'th element of Jv arising from this orbit as bJ2 v i = 0. 
This reasoning shows that the vector v stays localized 
under application of J. 

Second, we have to show that Jv still meets the zero- 
sum condition, so that a vector remains localized under 
iterated application of the Jacobian. For this purpose 
let us consider the connections of an orbit p with the 
nodes of an orbit q under the action of the Jacobian. One 
can always find a reordering of the variables such that 
all terms representing the links from nodes in p towards 
nodes in q appear in a single block in J. In the following 
we refer to such blocks as orbit connections and denote 
them by C qp . 

We first show that the sum over the elements of C qp v 
belonging to nodes in q is zero. More explicitly, we can 
write this sum as ^ 4 (C qp v) i , where i runs over the nodes 
in q (over all elements of the vector C qp v). However, this 
is identical to first taking the sum over the elements of 
each column of C qp and then multiplying by the vec- 
tor v. The column sums of C qp correspond each to the 
sum over all links originating from one node in p and are 
therefore identical by symmetry. Let this sum be c. We 
can then write J2i (C qp v) i = ■ Vj = 0, where i runs 
over nodes in q and j runs over all nodes in p. In other 
words, this shows that the sum of the input that a given 
orbit receives from any other orbit under application of 
the Jacobian on a vector in V sums to 0. Consequently 
the zero-sums conditions remain fulfilled for Jv such that 



v stays localized also under iterative application of the 
Jacobian matrix. 

Because the two conditions, localization and zero- 
sums, are met by Jv, V is an independent subspace under 
J, called the redundant subspace of the symmetry. We 
note that the redundant subspace of a symmetry of o or- 
bits and of size s has dimension o(s — 1). Further, the 
compliment of the redundant subspace W = V 1 - is also 
an independent subspace. The orthogonality of W and 
V implies that any vector w in W is constant on nodes 
belonging to the same orbit. 

Intuitively speaking the redundant subspace V cor- 
responds to cigenmodes of the Jacobian that charac- 
terize perturbations affecting nodes in an orbit anti- 
symmetrically. Thus the effect on nodes outside the re- 
spective symmetric structure cancels, which explains the 
localization. The subspace W correspond to eigenmodes 
of the Jacobian that characterize perturbations affecting 
nodes in the orbit identically. Thus the effect on nodes 
outside the symmetric structure does not cancel and cou- 
ples the orbit to the rest of the system. 

So far we have shown that the presence of a symmetric 
structure can be linked directly to the redundant eigen- 
values in the spectrum of the Jacobian. These eigen- 
values and the corresponding eigenvectors correspond to 
local eigenmodes governing the approach or departure of 
the system to or from the steady state. Therefore sym- 
metric structures provide a rare example of mesoscale 
motifs that have exactly prescribed implications on the 
system- level dynamics. 

For computing the redundant eigenvalues, one needs 
to diagonalizc the Jacobian only in the subspace V. The 
difficulty involved thus depends on the complexity of the 
symmetric structure, but not on the complexity of the 
embedding system. 

One can now imagine a mapping which collapses all 
nodes in a given orbit onto each other such that each or- 
bit is replaced with a single substitute node. Intuitively, 
this mapping should remove the redundant eigenvalues 
from the spectrum of J while leaving all other eigenval- 
ues unchanged. 

For confirmation we establish J for vectors inside W 
in terms of new variables w' p , each of which replaces all 
the (identical) variables of an orbit p. We recall that this 
substitution is possible because vectors in W are constant 
on each orbit. 

More explicitly, the elements of the vector Jw cor- 
responding to nodes of given orbit q are given by 
J2 P C qp lw' p , where p runs over all orbits, C qp denotes 
the orbit connection from p to q, and 1 is a vector that 
contains the entry 1 for each node in q. Because each of 
the row-sums of C qp represents the sum over the incom- 
ing links of one node in q from nodes in p, the entries of 
C qp l representing these row-sums are identical. We can 
therefore consistently replace the entries in C qp lw /p 
by w' q . Hence the connection from the variables w' p to 
w 'p [ s given by the row-sum of C qp . This defines the re- 
duced Jacobian matrix J' which contains at each position 
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this connection from w' q to w' p . Furthermore it defines 
the reduced (hypothetical) network of J' which contains 
one node for each replacement variable (w' p instead of 
the orbit p), confirming the intuitive assumption. 

These results can be summarized in a simple algorithm 
which reduces the network of the Jacobian to the network 
of the reduced Jacobian. Starting with the graph of the 
Jacobian J, we eliminate all nodes except one for each 
orbit. Links that pointed towards an eliminated node 
are removed and links that originated from an eliminated 
node are rewired to the remaining node from the same 
orbit. Finally, we simplify the network by replacing all 
links originating and ending in the same nodes by a single 
new link. The weight of this replacement link is the sum 
over the replaced ones. Repeating this procedure if new 
symmetries arise from the reduction we ultimately obtain 
the so-called quotient graph corresponding to the reduced 
Jacobian which contains no more symmetries. 

The spectrum of the reduced Jacobian matrix J' de- 
scribes the dynamics of the original system, except for 
the localized redundant part. However, if the redun- 
dant eigenvalues of J have negative real part, they do 
not change the dynamics qualitatively and analysing J' 
is equivalent to analysing J. 

Once the Jacobian has been equivalcntly reduced, one 
can try to construct a physical representation of the re- 
duced system. We illustrate this process for food webs in 
the next section. 



IV. APPLICATION TO FOOD WEBS 

The results in Sec. Mil imply that the nodes of each 
symmetry in a predator-prey network carry localized dy- 
namical modes inside the associated redundant subspace 
of the Jacobian matrix. One can imagine that reducing 
the food web model by collapsing the symmetric nodes 
removes the modes associated with the symmetry but 
leaves the remaining dynamical modes intact. Thus this 
reduction leads to an equivalent dynamical system, if the 
removed modes are stationary and stable, i.e. when the 
redundant eigenvalues have negative real parts. In this 
section we confirm this intuition by analysing different 
symmetric structures inside simple example food webs. 

We start by considering the so-called competitive- 
exclusion motif shown in Fig. [3^). In this motif two 
top-predators (populations 1,2) feed on the same prey 
(population 3). Using Ref. LL9| the stability of any given 
steady state is governed by a Jacobian of the form 



Vu Jl3 
J22 J-23 

.•hi J32 J33, 



(3) 



The diagonal terms reflect for each species the depen- 
dence on its own abundance of mortality, of primary pro- 
duction, and of the predatory losses and gains. The off- 
diagonal terms reflect the dependence of the predator- 
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FIG. 3. Examples for food webs that differ only by dynamical 
modes localized inside their mesoscale symmetries. Different 
colours denote different orbits. Two systems are marked as 
dynamically equivalent ( "=" ) if the spectra of their Jacobian 
matrices are identical up to the eigenmodes characterizing the 
mesoscale effects inside the symmetries. If nodes of a sym- 
metry have a common prey or predator with a node outside 
the symmetry, as e.g. in i), the generalized parameters of the 
common prey or predator in the reduced network are weighted 
which is shown in j) as a double-headed link. See text for more 
details. 



prey interactions on either the predator or the prey abun- 
dance. 

This motif contains a symmetry if the generalized 
parameters of population 1 and 2 are identical (such 
that J13—J23, J 11= J 22, J 31= J 32) and simple diagonaliza- 
tion yields the three eigenvalues A3 = Jn and Ai ; 2 = 
+ J33 ± VA), where A = (Jn - J33) 2 + 4J 3 iJ 3 i. 

Following the procedure outlined in the previous sec- 
tion, we can write the Jacobian of a smaller system 



J _ J "Al ^13 

\2Jsi J33 



(4) 



that exhibits qualitatively similar dynamics. Straightfor- 
ward diagonalization confirms that this reduced matrix 
J retains Eigenvalues A1.2, while the eigenvalue A3 cor- 
responding to the mode associated with the mesoscale 
symmetry has disappeared. 

We note that the reduced matrix J can be interpreted 
as the Jacobian of a predator-prey system. This is not 
trivial because the terms such as J13 and J31 are not 
independent, but can be expressed as functions of the 
same set of generalized parameters. Using the detailed 
representation of these terms [13, [HJ we verified (c.f. [5} 
that J corresponds to a predator-prey system, in which 
the prey is described by the same parameters as species 1 
in the original system and the predator by the parameters 
of species 3. It is remarkable that the factor 2 in J that 
arises from the reduction is indeed needed to make this 
direct correspondence possible. 

The dynamics of the reduced (predator-prey) system 
is equivalent to the original (competitive-exclusion) sys- 
tem if the dynamics in the direction of the redundant 
eigenvector is stable, i.e. if the eigenvalue A3 that dis- 
appears in the reduction has a negative real part. For 
the competitive-exclusion motif this is typically the case 
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when the predators suffer from super-linear mortality 

While we have so far focused on the three node web 
for illustration, we emphasize that the same reduction 
can be applied in food webs of arbitrary size. In particu- 
lar, using the results from the previous section, it is not 
necessary to compute the complete spectrum of the Ja- 
cobian to apply the reduction. For instance, consider the 
competitive-exclusion motif inside the complicated food 
web in Fig. [3Ji) for which writing all the Jacobian ma- 
trix eigenvalues is very complicated. Even in this case 
the results from Sec. IIIII guarantee that the the compet- 
itive exclusion motif that appears in the web carries the 
redundant eigenvector associated with the eigenvalue A3 
found before. Removing the associated dynamical modes 
leads to the Jacobian matrix of the reduced food web in 
Fig. [3J;), which is dynamically equivalent if A3 < 0. 

As a second example we consider the so-called 
apparent-competition motif shown in Fig. [3fc), i.e. one 
predator (population 1) feeding on two prey species (pop- 
ulations 2, 3). The Jacobian matrix is 

(J11 J12 Jl2\ 
J21 J22 J23 , (5) 
J21 J23 J22) 

where we have assumed that the prey species are de- 
scribed by identical ecological parameters. In contrast 
to the Jacobian of the competitive-exclusion motif this 
matrix contains the elements denoted by J23 that cou- 
ple the symmetric (prey) species. These so-called appar- 
ent competition terms arise because feeding on one prey 
species leads to a saturation of the predator that benefits 
both prey species. By contrast the two predators in the 
competitive-exclusion motif can generally be assumed to 
feed on the prey independently of each other. 

The dynamical mode localized on the two prey species 
is characterized by the eigenvalue A3 = J22 — J23 and 
the corresponding eigenvector (0, 1,— 1) T . Therefore if 
J22 < J23, which is generally the case in food webs, then 
removing the dynamical mode of the symmetry leads to 
a dynamically equivalent system. 

Following the procedure outlined in Sec. Illlj. we can 
write the Jacobian of the reduced system as 

t I Ju 2Ji2 \ , g n 

J = U j ™ + W ' (6) 

which retains the eigenvalues of J except for A3 . Compar- 
ing J to the Jacobian matrix of the system in Fig. [3b) we 
find that J describes a predator-prey system with iden- 
tical generalized parameters (c.f IB")) . 

The remaining examples in Fig. [3] can be treated along 
the same lines as above. Let us therefore comment only 
briefly on two new features that appear in these exam- 
ples. First, when aggregating symmetric species that 
share a common predator with other species that are not 
part of the same orbit then care has to be taken that 
the aggregate species enters with correct weight. For in- 
stance in Fig. if we assume that the top predator 



receives two third of its biomass from the two symmetric 
species. Then in the reduced system the species replac- 
ing the orbit has twice the weight in comparison to the 
non-symmetric prey on the right. 

Second, we note that symmetries consisting of more 
than one orbit may exhibit more complicated localized 
dynamics. For instance in the food web shown in Fig.[3Ji) 
the symmetric structure consists of two interconnected 
orbits. Diagonalization yields the redundant eigenvalues 
Al,2 = i(Jii + J 33 ±%/A) with A = (Jn-J 3 3) 2 + 4J 3 iJi3, 
where we used the indices 1 and 3 to denote a top- 
predator and its prey, respectively. The redundant eigen- 
values can thus form a complex-conjugate eigenvalue 
pair, potentially leading to oscillatory instabilities. 

We emphasize that symmetries present in the steady 
state do not generally persist when the system departs 
from the this state. Dynamics that is generated in a 
smymmetric mesoscale motif can therefore spread and 
impact the entire network. For instance in the example 
discussed above the oscillatory instability can be inferred 
solely from the analysis of the mesoscale motif, wheras 
the resulting oscillations may affect the entire network. 

In summary every mesoscale symmetry is associated 
with one or more dynamical modes that originate in the 
symmetric motif. If the dynamical modes are stable, then 
they can be removed from the system without changing 
the remaining dynamics by replacing each orbit of the 
symmetry by a single node. For this reduced system 
the Jacobian can be obtained according to the rules de- 
scribed in Sec. IIIII The reduced Jacobian can in general 
be interpreted as describing a smaller food web, however, 
potentially with altered parameter values. For the simple 
cases considered in Ref. [T3 the parameters of the origi- 
nal predator-prey network and those of the reduced net- 
work are so closely related that the dynamically equiva- 
lent systems could be spotted empirically. However, this 
is much harder for complicated topologies, as e.g. in 
Fig. [3jl), and for more elaborate models. Even in such 
complicated cases the aggregated Jacobian can be com- 
puted relatively straight forwardly with the procedure 
proposed here. Notably the complexity of this reduction 
increases with the size of the symmetric structure, but is 
independent of the size of the embedding network. 



V. CONCLUSION 

In this paper we showed that a previously observed 
equivalence in the dynamics of food webs [13] is a conse- 
quence of localized dynamical modes originating from on 
mesoscale symmetries. 

We extended a recent result [lO] to show that cer- 
tain dynamical modes localize on symmetric structures 
in food webs. In cases where such mesoscale structures 
are present these dynamical modes are independent from 
the network outside the symmetry. If the localized modes 
are stable, then aggregating the symmetric structures in a 
food web as described in Sec. lIIII removes them while leav- 
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ing the remaining dynamics intact, leading to a smaller 
but dynamically equivalent system. 

In principle the findings of this paper can be used to 
formulate reduced models for certain food webs. In prac- 
tice this application of the results will be of little impor- 
tance as removing symmetries leads, at best, to a rela- 
tively mild aggregation. 

However, from the point of view of physics, symmet- 
ric motifs in networks, provide an example of mesoscale 
structures that have distinct and exact dynamical im- 
plications for the dynamics of the network as a whole. 
These implications can be explored even in very large net- 
works because computation of the corresponding local- 
ized eigenvector and eigenvalue do not require the com- 
putation of the complete spectrum of the Jacobian. 



Appendix A: Reduction of the Jacobian matrix of 
the Competitive-Exclusion Motif 

In this appendix we show the explicit reduction of the 
Jacobian matrix for the competitive exclusion motif in 
Fig. [3^i) in terms of generalized modeling parameters of 
the model presented in Ref. [l^ and l32l Furthermore we 
show that the reduced Jacobian matrix J is indeed iden- 
tical to the Jacobian matrix of the simple predator-prey 
system in Fig. [3)3). In terms of generalized modeling pa- 
rameters ([3]) becomes 

(a x 0\ M-Mi 7i \ 

J= a x 7i , 

\0 a 3 J \ $3-71/ 

where is the biomass turnover rate of species i, ^ x is 
the sensititivity of predation to the abundance of both 
predators, fix is the sensitivity of mortality to predator 
abundance of both predators, $ 3 is the sensitivity of pri- 
mary production to the abundance of prey and 71 is the 
sensitivity of predation for both predators to the abun- 
dance of prey. The factors of i in the last row arise 
because each of the species 1 and 2 is responsible for half 



of the biomass loss of species 3 also their effect on species 
3 is weighed by this importance. 

Using the results in Sec. IIII1 the reduced Jacobian ma- 
trix (c.f. 0])) is 




ai 0\ / $1-/11 71 \ 
a 3 J \2 $ 3 -7iy ' 



This is identical to the Jacobian matrix for a predator 
prey-system in which the predator is described by the 
generalized parameters ax, $1, 71, fix and the prey species 
by $ 3 - 



Appendix B: Reduction of the Jacobian matrix of 
the Apparent-Competition Motif 

In terms of generalized modeling parameters, the Ja- 
cobian matrix in ([5]) of the apparent-competition motif 
in Fig. [3fc) is 

fax 0\ /*i-Mi 571 571 \ 

J= a 2 -$1 -V . 

\Q a 2 ) \ _2s^f $2-^/ 

where the parameters are the same as in[Al The factors of 
\ in the matrix arise because each of the two prey species 
is responsible for half of the biomass gain of species 1. In 
comparison to the Jacobian of the competitive-exclusion 
motif, the matrix contains the elements — 71 ~ 1 that cou- 
ple the symmetric (prey) species. These are the apparent 
competition terms arising because feeding on one prey 
species leads to a saturation of the predator that bene- 
fits both prey species. 

Using the results in Sec. lIIIl the reduced Jacobian ma- 
trix is 

- fax 0\ f^x-tii 2± 7 i \ 
V a 2 J \ -*i $ 2 -7iy ' 

Therefore J is also identical to the Jacobian matrix of the 
simple predator-prey system in which the predator is also 
described by the generalized parameters ax, $1,71,^1 
and the prey species by $3. 
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